Pressure induced orientational glass phase in molecular para- hydrogen. 
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We propose a theoretical description of possible orientational glass transition in solid molecular 
para-hydrogen and ortho-deuterium under pressure supposing that they are mixtures of J = and 
J = 2 states of molecules. The theory uses the basic concepts and methods of standard spin-glass 
theory. We expect our orientational glass to correspond the 11' phase of the high pressure hydrogen 
phase diagram. 

PACS numbers: 



I. INTRODUCTION 

The high pressure phase diagram of solid hydrogen and 
its isotopes is a fascinating subject of investigation during 
recent decades. Although the p — T phase diagrams of 
ortho^ H2 and para— D2 under not too high pressures are 
well understood, the same can not be said about para — 
H2 and ortho — D2 and about extremely high pressures 
[see, e.g., the reviews RefslH, 0, 

At low temperature and ambient pressure pure solid 
0—H2 and p~ D2 consisting of molecules with orbital an- 
gular moment J = 1 crystallize in hep lattice with rotat- 
ing molecules on the lattice sites. At lower temperature a 
transition takes place to the phase with the orientational 
long range order (LRO) of antiferroquadrupolar type. 
This transition is accompanied by a structural transi- 
tion to fee lattice. The transition temperature increases 
with pressure depending on intcrmolccular distance R as 
so indicating that it is due to the anisotropic electric 
quadrupole-quadrupole interaction (EQQ). This picture 
remains valid up to c 0.55, where c is the concen- 
tration of moment bearing molecules. At intermediate 
concentrations (c > 0.12) NMR experiments have been 
interpreted in terms of the freezing of the orientational 
degrees of freedom and the transition to the quadrupolar 
glass phase. [l|, 0, 01 At the lower c concentration there is 
no orientational ordering. 

Molecular p — H2 and o — D2 also crystallize in hep 
structure. However, at low pressure they remain in this 
structure up to OK. This phase has no orientational or- 
der (J = 0) and it is called phase I (or LP phase). At 
higher pressures (^llOGPa in p — H2 and ~28GPa in 
o — D2) solids transform to orientationally ordered bro- 
ken symmetry phases (phase II or BSP). [H, 0, @|- The 
possibility of orientational order in systems of initially 
spherically symmetric molecule states is due to the in- 
volving of higher order orbital moments J — 2, 4... in the 
physics under pressure. The crystal field of the neighbors 
perturbs the molecular wave functions and one can gain 
in overall energy if the anisotropic EQQ interaction be- 
tween nonzero quadrupolc moments is included. One has 
to keep in mind that EQQ interaction rapidly increases 
with increasing pressure. The long range orientational or- 
der appears abruptly at a fixed value of pressure through 



the first order phase transition just as it takes place in 
ortho-para mixtures when the concentration of moment 
bearing molecules achieves certain fixed value. 

Goncharov et al. investigated the high-resolution 
Raman spectra of almost pure o — D2. The authors in- 
dicate that in the intermediate pressure range between 
the phases I and II the ordering is incomplete and orien- 
tational frustration takes place. They further speculate 
that this intermediate IF phase exhibits glassy behavior. 
Phase IF persists for a narrow pressure range (~ 2Gpa) 
and has abrupt boundaries. 

It seems obvious that the T—p phase diagram contain- 
ing I, IF and II phases can be considered in close analogy 
to the T — c phase diagram of ortho-para mixtures. For 
simplicity we imagine that p maps C2 (p) - the concentra- 
tion of molecules with J = 2, although one should take 
into account other anisotropic interactions to understand 
the results of the precise experiments [e.g., Raman scat- 
tering, RefQ- The first attempts to describe in such 
a way the long range order and the orientational glass 
phase in o — D2 and n — H2 on a microscopic theory level 
were done in Refsll, @, HO. 



II. J = 2 QUADRUPOLE GLASS MODEL. 

The purpose of this paper is to give a theoretical de- 
scription of the possible orientational glass transition in 
solid molecular para-hydrogen and ortho-deuterium un- 
der pressure supposing that they are mixtures of J = 
and J = 2 states of molecules. The theory uses the ba- 
sic concepts and methods of standard spin-glass theory. 
We expect our orientational glass to correspond to the 
IF phase of the high pressure hydrogen phase diagram. 

It is well known that the number ofj = 0-^ J = 2 
transitions increases rapidly with increasing pressure [see, 
e.g. , Refsd, i, i, S 0, \m EJ, E [3 ■ The anisotropic in- 
teraction potential and the crystal field grow rapidly with 
increasing density. The energy of the many-body system 
can be lowered by taking advantage of the anisotropic in- 
teractions. So, the single molecule wave functions are no 
longer spherical symmetric but are rather admixtures in- 
cluding higher order excitations. This admixture is prob- 
ably responsible for the decrease of the critical conccn- 
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tration of J = 1 molecules for the LRO transition in 
ortho-para samples at high pressure the deficit of mo- 
mentum bearing molecules being compensated by J = 2 
molecules. It seems that the weak dependence of the 
transition to the III (or A) phase on the ortho-para com- 
position can be described in an analogous way. 

The rough estimation ofj = 0— >J = 2 transition 
probability ^ can be done using quantum mechanical 
perturbation theory consideringthe field of the nearest 
neighbours as the perturbation: [8| 



/25 r \ fRo^'' 



(1) 



where B is the rotational constant, R is the intermolecu- 
lar distance, T/B = 0.011 in H2 and 0.028 in 1)2. d The 
factor 3.8 has a geometrical nature and it corresponds 
to hep lattice. Using the compressibility data we obtain 
the pressure dependence of ^. This dependence is very 
strong. If we attribute to the probability ^ the mean- 
ing of the concentration of moment bearing molecules we 
see that the position and the width of IF phase qualita- 
tively coincide with that of quadrupolar glass in ortho- 
para mixtures. For example, ^ = 0.1 at 40 Gpa for o— D2. 
The isotope dependence is also roughly correct. 

The admixture of J = 2 states causes the frustration 
in the case of hep lattice. As to disorder, it is not ob- 
vious that one can think of J = 2 impurities as of the 
quenched disorder even at low temperature. Neverthe- 
less, it seems possible to consider the whole ensemble 
of physically achievable realizations of mixed states as a 
convenient background [see, e.g., RefflU for the formu- 
lation of an orientational glass model in the spirit of the 
spin-glass theory, [l^ 

Here we present two theoretical models of possible 
quadrupolar glass with J = 2. The first one is a general- 
ization of well known Sherrington-Kirkpatrick[l7j spin- 
glass and it is analogous to the model[ll| which describes 
well the quadrupolar glass in ortho-para mixtures. 0, Q 
The second model is a generalization of the so-called "p- 
spin glass" and it is probably more adequate for high 
pressures when many particle interactions can play an 
important role. We consider the case p = 3 in detail. 
The results obtained for the models differ: in the three- 
site model the discontinuities in the specific heat and in 
the glass order parameter as the functions of the temper- 
ature do appear. We hope that future experiments will 
discriminate between these models. The essential feature 
of the obtained intermediate phase in both models is the 
coexistence of the orientational glass with the long range 
orientational order as it is seen in the experiment. Q 



III. QUADRUPOLE GLASS WITH 
TWO-PARTICLE INTERACTION 

As the first model of the quadrupole glass we will con- 
sider a system of particles on lattice sites i, j with random 



truncated EQQ Hamiltonian 

H — — ^ ^ JijQiQj 



(2) 



Here Jij are random interactions distributed with the 



Gaussian probability 



N 



P{Jij) = —?==^ exp 



27rJ 



2J2 



(3) 



where the factor N insures the sensible thermodynamic 
limit. 

Operator Q ~ [3 J^^ — J( J + 1)] is the axial 
quadrupole moment of the hydro gen molecule in the 
space J = const; TrQ = 0. In RefllSl quadrupole glass 
freezing in ortho-para mixtures has been considered on 
the base of the Hamiltonian ^ in the subspace J = 1 
with — 0, ±1 and Q = 3 J^^ — 2, so that 



Qli) — 2 — Q(i 



(4) 



Now the Hamiltonian ^ will be considered in the sub- 
space J = 2 with Jz = 0, ±1, ±2 and 



Q = l [3J.' - 6] 



(5) 



Let us emphasize that the model with J ~ 2 differs es- 
sentially from that with J = 1 because of the different 
operator algebras. For example, now the operators Q are 
5x5 diagonal matrices and instead of (|4]) we have 



q3 = 4 + 4g - 



(6) 



Following the standard methods of the spin-glass the- 
ory [see, e.g., and using the replica technique, 
we can express the disorder-averaged free energy of the 
system in the form [see also Refs 
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-Y.iq'^'?-\nTr[Q.yexpe\, (7) 



a>0 



where 



(8) 



Here t = J /T and the numbers a and (3 label replicas. 

The standard conditions for the free energy saddle 
point determine the glass order parameter q and the reg- 
ular order parameter x [the average quadrupole moment] , 



(9) 
(10) 
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FIG. 1: (Color online). The temperature dependence of 
the order parameters for the quadrupole glass with the two- 
particle interaction. The replica symmetry breaking occurs 
at the temperature corresponding to A(RS)ropi = 0. 
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FIG. 2: The heat capacity of the quadrupole glass with two- 
particle interactions as a function of temperature. The arrow 
marks the temperature corresponding to the boundary be- 
tween RS and more stable IRSB solution. 



and also the auxiliary order parameter 

= ((Q")'>e, 

where 



Tr 


"(...)exp (e) 


Tr 


cxp 









(11) 



(12) 



In the repUca symmetric (RS) approximation, when 
all qap arc equal, the expression ^ for the free energy 
becomes 



Frs = -NT <; - + InTr ( cx^^Ors 



Here 



(...) 



dz 



zty/qQ + t 
(. . .)exp 



dz^ {...). 



(13) 

(14) 
(15) 



Using the extremum conditions for the free energy, 
Eq. ()13p . we obtain the equations for the order param- 
eters: 



1 - (Q)L - 

P 



Bus- 



(16) 
(17) 
(18) 



The results of numerical solution of the Eas. (|16llT5|) are 
illustrated in FiglT] One can see that there is the coex- 
istence of glass and LRO. The RS glass order parameter 



and the RS average quadrupole moment grow continu- 
ously on cooling and are nonzero even at arbitrary high 
temperature. The absence of the zero solution of Eqs. pGI 
[TH]) follows from the fact that TrQ^/^^+i =^ 0, A: 1, 2, . . . 
[see Refl20l for details] . The orientational contribution to 
the heat capacity 



N 



d{l/t) 



(19) 



has a broad maximum at the temperature slightly lower 
than that of the instability of the RS solution. All these 
features are common for the quadrupole glass with J = 2 
and J 1 [see Refsi, 

The replica symmetric solution is stable unless the 
replicon mode energy A is nonzero. [2 1| For our model 
we have: 



'^(RS)rcpl — 1 — i 



(20) 



At the temperature Tq defined by the condition 
-^(RS)rcpi = the RS solution becomes unstable and one 
needs to look for the solutions with the broken replica 
symmetry ( RSB ) . Using the standard procedure [see, 
e.g., RefUJ, we perform the first stage of the replica 
symmetry breaking (IRSB) according to Parisi [n repli- 
cas are divided into n/m groups with m replicas in each] 
and obtain the free energy in the form [with q""^ = ri if 
a and /? are from the different groups and = ri + v 
if a and (3 belong to the same group] 



IRSB 



NTx 

rr ri + w - p- 
m— + (1 — m) 



(21) 



dz^ln / ds^ 



Tie 



SlRSE 
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where 

eiRSB zt^Q + stVT>Q + t' ^ ~ ~ "" g^ (22) 

The extrenium conditions for Fi^gg yield the equations 
for the glass order parameters ri and v, the regular order 
parameter x, the additional order parameter p, and the 
parameter m [see Appendix |^ where Ou^sb is given by 
Eq.(l22l) for Z = 2]. 

To estimate the form of the IRSB solution near the bi- 
furcation point, To, at which it ceases to coincide with the 
RS solution [i.e., in the neighborhood of Tq], we expand 
the expression for the free energy ([7]) up to the third or- 
der inclusively, assuming that the deviations Sq"^ from 
qjis are small. In fact we expand the argument of the 
exponent: 



a>l3 



ypE(^")'+*'tosE^"Q^ (23) 

with t = to + At. Using the formulas of Appendix [BJ we 
obtain 



AF 

NT 



{1 - t^W) [r - (m - l)vf + 



{m-l)} + jL[r - {m-l)vY 



(24) 



[r - {m - l)vf + 

D [r — (m — l)v] v'^m{m — 1) — 

Bsv^m^im - 1) + BiV^m{m - 1)(2to - 1)} + ... 

where t = to + At, r = ri — qRs^ ^md the expressions for 
the parameters W, L,C, D, B^, are given in Appendix 

Using the extremum conditions for the free energy 
and the fact that L\t=to ^ 0, we obtain the branching 
condition r — (m ~ l)f = + o(At)^, i.e., the condition 
that there is no linear term in the glass order parameters. 
There is no other linear term because (1 — t^VF)|t=tQ = 
'^(RS)rcpi|t=to = at the bifurcation point. Finally, we 
obtain: 



to t% dW , 



\t = to 



At = tl [-Bi + m{-B3 + 2Bi)] v, 



{2m- I] 



to 
2 



tJdW. 



\t = to 



At 



t^U2m-l) -B^ + m{-B3 + 2Bi) 



mim-l){-B3 + 2Bi) \v, 



(25) 



where B^ and B4 are taken at T = To. So, 

At; r= (to - l)v, (26) 



in the neighborhood of Tq, where the IRSB solution ap- 
pears and 



B4 
B3 



(27) 



at the branch point Tq. 

Let us notice that all the obtained expressions hold for 
Hamiltonian Q , where Q is the arbitrary diagonal oper- 
ator such that TrQ = 0, TrQ'^ and T|t=to are nonzero. 

For our model with Q defined by Eq.®, it follows 
from above formulas that to = 0.25 and there is no jump 
in the order parameters at the point where the IRSB 
solution appears [as usually when m < 1]. The solutions 
of the equations giving the extremum conditions of (|2ip 
are presented in Fig[T] The orientational order and the 
glass regime coexist and grow smoothly on cooling even 
through RS - IRSB transition. In addition, the curve for 
the heat capacity changes a little in passing from the RS 
to IRSB solution (FigH]). 

The IRSB solution is stable above the temperature 
determined by the solution of the condition, A(iRgB)ropi = 
[see FigU, 



-^(IRS)rcpl — 1 — t X 



dz 



G 



Ids 



G 



Tr e^'i^^s 



(<3^)einsB - ((Q)eiRSB)^ 



JdsO 



(28) 



At T < To we have the nonergodic state. At T ^ one 
expects the full replica symmetry breaking (FRSB). 



IV. p-SPIN GLASS LIKE QUADRUPOLE 
MODEL. 



In this Section we consider a generalized p-spin inter- 
action spin glass model — p-quadrupole model. The 3- 
quadrupole case will be considered in detail (see below). 
The model of Sec lIIIl is 2-quadrupole model. The p-spin 
glass model with the random interaction of p Ising spins 
was considered in a large number of papers [see, e.g., 
RefsH^, m, m, [2^, |2^ and serves as a generic model for 
investigation of glasses without reflection symmetry. 

Now we consider I- quadrupole model described by the 
Hamiltonian 



H 



^ ^ Jii...ii QiiQi 



(29) 



ii<i2...<ii 



where i = 1,2, ...N,and Q is defined in Sec lIIII The cou- 
pling strengths are independent random variables with a 
Gaussian distribution 



P{Jn...ii) = rn-~ exp 



IIttJ 



(30) 
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FIG. 3: (Color online) Order parameters evolution with the 
temperature for 3-quadrupole model. The transition RS - 
IRSB takes place at the point defined by the condition m = 1. 
Glass order parameter v has a jump at this point. 



Using the replica approach we can write the free energy 
averaged over the disorder [see for details Reff2^ in the 
form: 

{F)jINT ^ hm - max | ^(K)' + ^ ^"(^^")+ 

a Q 

-T E + E ^"''9"'' - l^Tr{Q»} exp4 . (31) 



where 



Q>/3 a 



(32) 



The extremum in Eq. pip is taken over the physical 
order parameters and over the corresponding Lagrange 
multipliers, A"'^ and /.t". So the saddle point conditions 
give the glass order parameter 



(33) 



the regular order parameter [average quadrupole mo- 
ment! 



the auxiliary order parameter 
and the parameters 



(34) 
(35) 



X'-P =^-l{q"Pf~^\ = ^;(p")('-i). (36) 




FIG. 4: The heat capacity of the quadrupole glass with 3- 
particle interactions as a function of temperature. There is a 
jump at the RS - IRSB transition. 



Using the standard procedure [see, e.g., ReflTol we can 
obtain RS and IRSB expressions for the free energy. Let 
us write Fij^sB for the case I = 3. The free energy Fas 
can be obtained if one put v = 0. 

FiRSB = -NTx 



2 ^ ; r, 



(37) 



1 

771 



dz^ln / ds 



„G 



Trcxp (^OiRSB^ 



Here 



^IRSB = Zt\j Ci- 



st \ 



'3[(ri+i;)'-ri2] 



Q- 



,.3[p2-(n+.)1^.^ (38) 



The extremum conditions for FmsB yield the equa- 
tions for the glass order parameters, ri and v, the regular 
order parameter x, the additional order parameter p, and 
the parameter 777 [see Apendix |^ where Oirsb is given 
by Eq.® for I = 3]. 

It is easy to show that the corresponding condition 
•^(RS)rcpi — docs not determine a physical solution in 
the vicinity of the bifurcation point Tq: one passes to the 
unphysical free energy branch. In fact, the transition RS 
- IRSB takes place at the point Ti > Tq defined by the 
marginality conditions. At this point m = 1 and Frs = 
FiRSB- There are no discontinuities in free energy. The 
order parameters r, x and p are continuous and v has a 
jump [sec Fig[3]. When the temperature is decreased 777 
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becomes smaller than one and the IRSB solution leads 
to a larger (preferable) free energy than the RS solution. 

The corresponding discontinuity occurs also in the heat 
capacity, 



(IRSB) 



mrf + (1 — m) (ri + v)^ ^ 



(39) 



The form of the curve for the heat capacity, FigHl is 
analogous to obtained in RefH^ for the spherical p-spin 
model. 

The IRSB solution can be stable above the tempera- 
ture T ~ T2 determined by the second solution of the 
condition A(ip{gB)rcpi = 0, see FigI3J 



A 



(IRSB)rcpl 



dz 



G' 





Tr e^"i«sB 


m 


l(Q^)6Iirsb 


- ((<3)eiRSB)^} 


JdsO 




m 



(40) 



At the point T2 a transition to FRSB-statc or to a stable 
2RSB-state may take place. 



glass order parameter as functions of the temperature do 
appear. We hope that future experiments will discrimi- 
nate between these models. The essential feature of the 
obtained intermediate phase in both models is the co- 
existence of the orientational glass with the long range 
oricntational order as it is seen experimentally.! 
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APPENDIX A: 

The equations for IRSB glass order parameters ri and 
V, the regular order parameter x, the additional order 
parameter p, and the parameter m. 



V. CONCLUSIONS ?i 

In this paper we give a theoretical description of possi- 
ble orientational glass transition in solid molecular para- 
hydrogcn and ortho-deuterium under pressure supposing 
that they are mixtures of J = and J ~ 2 states of 
molecules. The theory uses the basic concepts and meth- 
ods of the standard spin-glass theory. We expect that ^ + ''1 
our orientational glass corresponds to the IF phase of 
the high pressure hydrogen phase diagram. 

We present two theoretical models of possible 
quadrupole glass with J = 2. The first one is a gener- ^ = 

alization of thewcU known Sherrington-Kirkpatrick spin- 
glass. The second model is a generalization of so-called 
"p-spin glass" and it is probably more adequate for high 
pressures when many particle interactions can play an 1 
important role. We consider in detail the case p = 3. 
The results obtained for two models differ: in three-site 
model the discontinuities in the heat capacity and in the and 
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Tr e'^iRSB 
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Tr e^iHSB 
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(A3) 



, (A4) 



m-(/-l)[(rl+^0 -(ri)'] = 



In 



Tr e' 



BlRSE 





Tr e^"i«sB 


In 


Tr e^'i^ss 




Tr 


Q01RSB 
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APPENDIX B: 

The only nonzero sums are 



and 



lim-V' Sq'^^Sq^'^Sq^" ^ {m-l){m^2W ~2,m{m-l)iji^ + 2n?i^-, (B2) 
lim i V' = (m- 1)27^3 _ j^Wj^^2 ^^^2^N ^^2^3. 

limly' V* = {m-lfri'^ + ?,TT?{m-l)ri^'^- imim - if if - nV^ , (B4) 

where 77 = r + w and ^ = r. The prhiic on the sum means that only the superscripts belonging to the same 6q are 
necessarily different in 



APPENDIX C: 



The formulas used to calculate the parameters are 



W = {QlQl) - 2{QlQ2Q3) + (Q1Q2Q3Q4); (CI) 

L = ~{QlQ2Q3) + {QiQ2Q3Qi); (C2) 

C = ~{B2 + B'^) + 2B3 + B'^- B^- (C3) 

D = -3^3-5^ + 354; (C4) 



where 



B'2 = ^{QiQ2Q3Q4Q5Q6) - ^{QlQ2Q3Q4Q5) + ^{QlQ2Q3Q4) 



B3 = liQlQlQl) ~l{QiQlQ3Qi) -l{QiQ2Q3Q4Q5Q6) + 1^^^^ 

B'3 = -(Q1Q2Q3Q4Q5Q6) + \{QIQ2Q3QaQ5) - \{QIQ2Q3Qa) - \{QlQlQ3Qi) + \{QIqIQ3); 
Bi = ^{QiQ2Q3QiQ5Q6) - {QlQ2Q3QiQ5) + l{QlQ2Q3Qi) + l{QlQlQ3Q4) - ^(QjQaQs) + j^iQl'^'' 



2I1 



and 



Tr 



(g^)'^-(Q«)"...expS 



Tr [exp ^] 
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